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Abstract 

The charm and beauty quark production in deep inelastic scattering at low values of the 
Bjorken variable x is considered in the DGLAP improved saturation model. After fitting 
parameters of the model to the structure function F 2 , the heavy quark contributions F£ c and 
F$ are predicted. A good description of the data is found. Predictions for the longitudinal 
structure function Fl and the diffractive structure function F® are also presented. 
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1 Introduction 



Among the color dipole models of deep inelastic scattering at small values of the Bjorken variable 
x, the saturation model [1,2] (GBW model) turned out to be especially successful. It was able 
to describe both the low x structure function F2 and the diffractive structure function F 2 D 
measured at HERA. From the theoretical side, the model has attracted a lot of attention since 
it grasps essential elements of parton saturation [3-5] incorporated in a relatively simple way. 
In particular, geometric scaling in DIS was predicted as a consequence of the assumptions on 
the saturation scale and scaling behaviour of the dipole cross section (which is a crucial element 
of the dipole models) [6]. With the advent of more precise HERA data [7-9], the GBW model 
needed improvement in order to provide better description of F2 at large values of the photon 
vir tuality (Q 2 > 20 GeV 2 ). This was done by Bartels, Golec-Biernat and Kowalski (BGK) 
by incorporating into the saturation model a proper gluon density evolving according to the 
DGLAP evolution equations [10]. This modification improves the small-r part of the dipole 
cross section. 

In the GBW saturation model the heavy quark contribution to Fi was considered in the 
form of the cc pair production. This element, however, is missing in the BGK improvement. 
The recent data from HERA [11-13] shows that the heavy quark contribution is of the size up 
to 30% and can by no means be neglected. Thus, the main goal of this analysis is to take into 
account heavy quark production in the DGLAP improved saturation model and confront it with 
the recent data. This analysis does not introduce new parameters to those already present. 
Once the parameters of the dipole cross section are determined from a fit to the total structure 
function F2, they can be used to predict the charm and beauty contributions F% c and F^, which 
are not fitted separately. In addition, the longitudinal structure function Fl and the diffractive 
structure function F® can also be predicted. 

It should be stressed that other dipole models are also successful in the description of the 
DIS data. These include the model of Forshaw and Shaw [14-16] based on the Regge-like ideas 
combined with the idea of saturation, the model of Iancu, Itakura and Munier [17] closely related 
to the theory of the color glass condensate [18-25] and the model of McDermott, Frankfurt, 
Guzey and Strikman [26] which is close in spirit to our analysis. Quite recently, a new analysis 
of Kowalski, Motyka and Watt [27] appeared which addresses a very important problem of the 
impact parameter dependence of the dipole scattering amplitude. For more details on these 
models see recent review [28]. 

The outline of this paper is as follows. In Section |2] we recall the main features of the BGK 
saturation model. In Section we describe fits with heavy flavours to the data on F2 from 
HERA. In Section we discuss main features of parton saturation, namely the critical line 
and saturation scale. Based on the results of the fits, in Section El we predict the charm and 
beauty contributions to F2 and the longitudinal structure function. In Sectional we discuss the 
diffractive structure function. Finally, conclusions are presented in Section 

2 The BGK saturation model 

In the dipole picture of DIS at small values of x [29-33], the virtual photon-proton scattering 
cross sections Oj. £ from transverse (T) and longitudinal (L) photons are given as a convolution 
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of the photon wave functions i^t,l an d the dipole cross section a, 



Nf r 

a^l(x,Q 2 ) = X) / d 2 rdz\if) T ,L( r > z iQ 2 'i m f> e f)\ 2 &{x,r) , (1) 
f=i J 

where / denotes flavour of the quark-antiquark pairs (dipoles) interacting with the proton (m/ 
and ef are quark mass and electric charge, respectively). The integration is performed over the 
qq transverse separation vector r and the light-cone momentum fraction of the virtual photon z 
carried by the quark or antiquark. The proton structure function F2 is proportional to the sum 
of the two cross sections 

^ = 7^-(4 P + ^ P )- (2) 



In the BGK model [10] the dipole cross section was assumed in the form 

a{x, r) = cr I 1 - exp I — I > , (3) 

which recovers the well known result of perturbative QCD [34] that for small dipole sizes the 
dipole cross section is proportional to the gluon distribution g(x,/i 2 ) with the scale /j? ~ 1/r 2 . 
Indeed, expanding j3J) in powers of r 2 , we find to the first order 

7T 2 

a(x,r) ~ —r 2 a s xg(x,n 2 ). (4) 

The feature that <t(x, r) vanishes in the limit r — > is called colour transparency and results 
from gauge invariance of QCD. The gluon distribution evolves with the scale according to the 
leading order DGLAP evolution equations. Since the model is applied in the low x region, 
quarks are neglected in the evolution. Inspired by the MRST parametrisation [35], the initial 
gluon distribution for the DGLAP evolution is taken in the form with two parameters, A g and 
X g , to be determined from fits to data 

xg(x,Q 2 Q ) = A g x x "(l-xf 6 . (5) 

The initial scale Qq is conventionally set to 1 GeV 2 

The behaviour of the dipole cross section (j3J) at large r has not been changed with respect 
to the GBW model in which the dipole cross section saturates at the maximal value do (also 
to be determined from fits to data). This key assumption is in agreement with unitarity of the 
cross section cr 7 p . The large dipole sizes, however, are the cause of troubles for the scale of the 
gluon distribution in Eq. © since fi may become smaller than Aqcd- In order to avoid this we 
freeze the scale at some minimal value ^0 3> Aqcd- Thus we assume that 

fj 2 = ^ + ^l, (6) 

which introduces two additional parameters, C and Hq, to be found from fits to data. For r — > 
the presence of /j,q becomes irrelevant and we recover the relation fj? ~ 1/r 2 . 

On the whole, the BGK model has five parameters to be fitted: the dipole cross section 
bound <To and the four parameters of the gluon distribution: A g , \ g , C and /jLq. In comparison 
to the GBW saturation model, there are two new parameters introduced by the scale 
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3 Fit description 



We performed fits with the charm and beauty contributions in the sum in Eq. (^Q) using recent 
data on the proton structure function Fi from HI [7] and ZEUS [8,9]. We considered only the 
data points with x < 0.01 and Q 2 > 0.04 GeV 2 . The number of experimental points in such a 
case was equal to 288. The statistical and systematic errors were added in quadrature. Moreover, 
the HI data were multiplied by a factor 1.05 to account for slightly different normalisation of 
the HI and ZEUS data sets. Similarly to the analyses [1,10], we also modified the argument in 
the dipole cross section a(x,r) in the heavy flavour contributions, 

/ 4m 2 , \ Q 2 + 4m 2 f 

x ^ x \ l + w) = WTw^^ (7) 

where W is invariant energy of the "f*p system. This modification takes into account that the 
heavy quark masses play the role of the hard scale when Q 2 <C m 2 b . 

We took typical values of the heavy quark masses, m c = 1.3 GeV and = 5.0 GeV, but 
the light quarks stayed massless. Therefore, we excluded the photoproduction point Q 2 = 
from considerations since in the dipole models <r 7P depends logarithmically on the quark mass 
in the limit Q 2 — ► 0. A comment is in order at this point. We performed preliminary fits with 
the light quark mass m q = 140 MeV taken from [1]. However, we did not find the value of the 
photoproduction cross section measured at HERA (174 fib) [36]. Thus the description of pho- 
toproduction would require fine-tuning of the light quark mass in the region where perturbative 
QCD, which is the basis of the dipole models, is not valid. Nevertheless, in the case of the heavy 
flavour production, when the quark mass is the hard scale, the predictions for Q 2 = could be 
done. 

Before we present the results of our fits, let us mention that the GBW model with heavy 
flavours fitted to the new data [7-9] gives x 2 / n df = 2.32. As it has already been pointed out 
in [10], this rather poor fit quality results from the lack of the proper DGLAP evolution of the 
gluon distribution in the GBW model. 

We performed two fits with the dipole cross section ([3]). taking into account the charm 
and beauty contribution in addition to the three light quarks. In the first fit only charm was 
considered while in the second one both heavy flavours were present. We set the number of 
active flavours in a s to 4 and 5, respectively and the value of Aqcd = 300 MeV in both cases. 
The fit results for the five parameters of the model, oo, A, X g , C and fj,Q, are presented in Tabled 
In the last row we recall from [10] the parameters of the BGK fit with the light quarks only. 

As we see, the value of x 2 / n df is still good for the fits with heavy flavours. The gluon 
parameters differ significantly from the light quark fit. In particular, the power X g is negative 
which means that the initial gluon distribution (JSJ) grows with decreasing x, in contrast to the 
fit with light quarks only when the gluon distribution is valence-like (\ g is positive). With the 
found gluon density, the total proton momentum fraction carried by gluons is around 25% at 
the inital scale Qq = 1 GeV 2 . 

In Fig. ^ we show the comparison of the dipole cross sections from the present analysis with 
heavy quarks (solid lines) and the BGK analysis [10] with light quarks only (dashed lines). The 
effect of heavy quarks is seen in the shift of the dipole cross section towards larger values of 
r, which means that for a given dipole size saturation occurs at lower x (higher energy). This 
effect was also observed in the GBW analysis [1]. 
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A, 


\ 


C 


Mo 


X 2 /ndf 


light + c + b 


22.7 


1.23 


- 0.080 


0.35 


1.60 


1.16 


light + c 


22.4 


1.35 


- 0.079 


0.38 


1.73 


1.06 


light 


23.8 


13.71 


0.41 


11.10 


0.52 


0.97 



Table 1: The parameters of our fits with heavy quarks to the F% data. The results of the BGK 
fit from [10] with massless light quarks only are given for the reference in the last row. 

It should also be mentioned that the presence of heavy quarks in the DGLAP improved 
model cures the pathological behaviour of the dipole cross section found in [10] in the fit with 
massless quarks (FIT 2). This is why we choose the BGK fit with massive light quarks (FIT 1) 
for the comparison in Fig. Q 



4 Critical line and saturation scale 

The shift of the dipole cross section towards larger values of r has direct impact on the position 
of the critical line which marks the transition to the saturation region in the (x, Q 2 )-plane. The 
line is defined by the condition 

&o(x,f)=aa , (8) 

with the constant a of the order of unity and the dipole size taken at its characteristic value 
f = 2/Q. In the GBW and BGK analyses the critical line was defined by the condition that 
the argument of the exponent in the dipole cross section equals 1, which corresponds to a = 
1 — e _1 ~ 0.63. With the form ©, the following implicit relation between x and Q 2 is found in 
such a case 

47T 2 

3(To q2 a s {^ 2 )xg{x,fi 2 ) = 1, (9) 

where the scale /i 2 = CQ 2 /4 + (Iq. This equation can be solved numerically to obtain the critical 
line shown in Fig. El as the solid line. The saturation effects are important to the left of this 
line. For the comparison, we also show the critical lines from the BGK and GBW analysis 
with light quarks only. We observe that the presence of heavy quarks shifts the critical line 
towards smaller values of Q 2 . This means that for a given Q 2 we need lower x in order to stay 
in the domain where the saturation effects are important. In other words, heavy quarks make 
saturation more difficult to observe at present and also future colliders, which is indicated in 
Fig. |21 by the acceptance regions of HERA and the LHC. 

An important element of the description of parton saturation in the GBW model [1] was the 
saturation scale, Q s (x) = Q s ox~ x , built in the dipole cross section 

&(x, r) = a {l- exp (- \r 2 Q 2 s {x)) } • (10) 

This form of the dipole cross section features the scaling behaviour , &(x,r) = a{rQ s {x)), which 
for massless quarks leads to geometric scaling of the total cross section, 

a^ = a^(Q 2 /Q 2 (x)), (11) 

found to a good accuracy in the small x HERA data [6,37]. Notice that the critical line in the 
GBW model is defined as a simple condition: Q 2 = Q 2 (x). 
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The BGK modification of the saturation model seems to abandon these important elements. 
Fortunately, this is not the case since in the fits with heavy flavours the value of C in the gluon 
scale © is small whereas /jLq ~ 1.6 GeV 2 is relatively large. It means that for not too small r, the 
dipole cross sections © effectively features the scaling at large values of r with the saturation 
scale proportional to the gluon distribution at the scale /i 2 , 

47T 2 

Q 2 s (x) ~ o — «s(Mo) x 9(x, Mo) • ( 12 ) 

Indeed, as we see in Fig. |21 which shows the dipole cross section as a function of the scaling 
variable r 2 Q 2 (x), geometric scaling is preserved for moderate values of r. It is broken, however, 
for small dipole sizes due to the DGLAP evolution of the gluon in the dipole cross section. 

It is important to notice that from the theoretical perspective, based on the theory of the 
color glass condensate (see [38] for most recent review), our considerations are confined to the 
mean field approximation which neglects the effects of fluctuations in the number of color dipoles 
in the proton wave function. These effects are important in the low density region of dipoles 
with very small sizes. For sufficiently small x (large rapidity Y = ln(l/x)), the fluctuations wash 
out geometric scaling leading to diffusive scaling [39-46]. The DIS cross sections in this case 
become functions of the variable ln[<5 2 /Q 2 (Y)]/ W DY rather than Q 2 /Q 2 (Y). The discussion of 
the new type of scaling is, however, beyond the scope of the presented analysis. 



5 Predictions for inclusive structure functions 

Once the parameters of the dipole cross sections (j3J) are determined from the fit to the F 2 
data, various inclusive structure functions can be predicted. In particular, the heavy quark 
contributions, F 2 C and F 2 b , can be found from the flavour decomposition in Eq. Q 

F 2 = F l 2 ight + Ff + Ff . (13) 

The dependence on flavour comes through electric charge ej and mass mj in the photon wave 
functions i/jt,l- For heavy quarks the mass also enters the dipole cross section through the 
modified value of the Bjorken variable ((TJ). This allows to compute the cc and bb photoproduction 
cross sections. For the HERA energy W = 209 GeV, we found 19.3 fib and 0.7 fib, respectively. 
To our astonishment, substituting the mass m q = 140 MeV for the three light quarks to the 
formula for cr 7P and performing then the photoproduction limit Q 2 — ► 0, we found 177 fib which 
agrees with the measured value 174 fib up to the experimental errors. 

The predictions for F 2 5 and F 2 b computed with the parameters from the first line of Table 1 
are presented as the solid lines in Figs. ^ and |S1 respectively. For the comparison, we present 
the predictions of the GBW model without the DGLAP modification (dashed lines). We see 
very good agreement with the data from HERA, both in the normalisation and the slope in 
x, in contrast to the GBW results which overshoot the data at large values of Q 2 . Thus, the 
presence of the DGLAP evolution in the BGK model is crucial for the correct predictions at 
large Q 2 . This can be understood by analysing the contribution of different dipole sizes to the 
heavy quark structure functions. For Q 2 3> 4m 2 , they are mostly sensitive to the small-r part 
of the dipole cross section (r ~ 2/Q) which is strongly modified by the DGLAP improvement. 

In Figure El we present the longitudinal structure function from our analysis (solid line) , 
plotted against Q 2 for W = 276 GeV. The experimental points represent the HI estimations of 
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Fl [7,47,48]. Reasonable agreement is found, however, the estimation errors are too large to 
draw firm conclusions. We also show the charm and beauty contribution to Fl (dashed line). 
We observe that in our analysis heavy quarks are important for large values of Q 2 while for 
Q 2 < 10 GeV 2 they may safely be neglected. 



6 DIS diffraction 

Deep inelastic diffraction is an important test of the dipole models [2,49,50]. The dipole cross 
section from the inclusive analysis can be applied to the description of the diffractive structure 
function. Following the approach presented in [2], the diffractive structure function is the sum 
of three contributions 

if (3) {xjp, 13, Q 2 ) = Ff + Ft + F™ , (14) 

which correspond to the diffractive systems composed of the qq and qqg dipoles with invariant 
mass M, produced by the transversely and longitudinally polarised virtual photons. Their 
interaction with the proton is described by the dipole cross section <t(xjp, r). Note that 

Q 2 + M 2 x 

Xp = WTW- 2 = -p (15) 

is substituted instead of the Bjorken gument of a. The longitudinal component F^ 9 

is suppressed by 1/Q 2 and is negligible in the region of small values of j3 = Q 2 /(Q 2 + M 2 ), i.e. 
for large diffractive masses M 2 ^> Q 2 , where the transverse qqg component dominates. 

The qqg diffractive amplitude was computed in [2] in the two-gluon exchange approximation 
with strong ordering of transverse momenta of the qq pair and the gluon. Each contribution is 
proportional to the inverse of the diffractive slope Bo for which we took the experimental value 
Bd = 6 GeV~ 2 . In addition, the qqg component is proportional to the strong coupling which 
we took a typical a s = 0.2. However, this is not quite clear which value of a s should be used, 
thus the normalisation of this term is rather uncertain. 

In Fig. [7| we show the comparison of the diffractive structure function F®^ computed in 
our analysis with the data from HI [51,52] and ZEUS [53,54]. However, these data are obtained 
for different upper limits on the mass of the proton dissociation system, My < 2.3 GeV for 
ZEUS and My < 1.6 GeV for HI. To account for this, the ZEUS data in Fig. is multiplied by 
the global factor 0.86 in accordance with [51]. At the same time we multiplied our predictions 
for the nondissociated proton by 1.23. This number, taken from [51], is the correction factor 
which takes into account the proton dissociation up to the mass My < 1.6 GeV. Thus, both the 
data and the theoretical predictions are normalised to the same experimental situation. The 
parameters from the first line of Table 1 were used in the dipole cross section for the computation 
ofif (3) . As we see from Fig. reasonable agreement with the data is found. 

In contrast to the inclusive case, the cc and bb contributions to the diffractive structure 
function are negligible due to the phase space restrictions for the diffractive system. The analysis 
of the diffractive ccg production, however, deserves more detailed considerations based on the 
seminal work [55,56]. 
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7 Conclusions 



We discussed the production of the charm and beauty flavours in the DGLAP improved satura- 
tion model [10]. We fitted the parameters of the model to the recent F% data from HERA, taking 
into account the heavy quark contributions in the theoretical formula for Fi- We found a good 
quality fit with x 2 /ndi close to unity. Thus we conclude that the successful description of the 
inclusive F2 data at low x which was found for the BGK model with light quarks is preserved 
when heavy flavours are considered. There are, however, some differences introduced by heavy 
quarks. First of all, the parameters vary significantly for the models with and without heavy 
flavours. This results in the shift of the dipole cross section towards larger values of the dipole 
sizes r with respect to the light quark case. As a consequence, the critical line in (x, Q 2 )-plane 
moves in the direction of smaller values of Q 2 which makes saturation more difficult to observe. 

The new predictions provided by our analysis concern the charm and beauty structure func- 
tions F^ and F^- We found very good agreement with HI and ZEUS data in all Q 2 bins. The 
significant improvement of the slope in x for high Q 2 with respect the GBW model is attributed 
to the DGLAP evolution. 

Other structure functions were also determined. For the longitudinal structure function Fl, 
we found agreement with the HI estimations. However, large estimation errors prevent from 
making more precise statement. The comparison with direct measurements, which are planned 
by HI and ZEUS collaborations, will be particularly interesting. For the diffractive structure 
function F% , we observe reasonable agreement with the HI and ZEUS data after a proper 
renormalisation of the data points. 

Finally, we discussed the issues related to the essential features of parton saturation like the 
saturation scale and geometric scaling. We showed that saturation scale is effectively present in 
the DGLAP improved model being related to the gluon distribution at the frozen value of the 
factorisation scale /Aq. 
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Figure 1: The dipole cross section in the BGK model with and without heavy quarks (solid and 
dashed lines, respectively) for x = 10~ 2 . . . 10~ 6 . 
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Figure 2: The dipole cross section in the BGK model with heavy quarks as a function of 
the scaling variable r 2 Ql(x) with the saturation scale given by Eq. (|12j) . Geometric scaling 
is preserved for moderate dipole sizes, while it is broken for small values of r because of the 
DGLAP modification. 
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Figure 3: The critical line in the (x, Q 2 )-plane from various saturation models indicating the 
position of the saturation region (to the left of these lines) . The shaded areas show the acceptance 
regions of HERA and the LHC. The latter region corresponds to the production of an object 
with the minimal mass squared Q 2 = 100 GeV 2 [57]. 
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Figure 4: Predictions for the charm structure function F^ in the BGK model with heavy quarks 
(solid lines). The predictions in the GBW model [1] are shown for reference (dashed lines). 
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Figure 5: Predictions for the beauty structure function i 7 ! in the BGK model with heavy quarks 
(solid lines). The predictions in the GBW model [1] are shown for reference (dashed lines). 
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Figure 6: The longitudinal structure function predicted in the BGK model with heavy quarks 
together with the HI estimations for various Q 2 at constant energy W = 276 GeV. 
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Figure 7: The diffractive structure function xpF D ^ (xp, (3,Q 2 ). The predictions of our analysis 
with heavy quarks are compared with the recent HI and ZEUS data. 
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